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ABSTRACT 

We report the result of the analysis of the light curve of a caustic-crossing binary-lens mi- 
crolensing event OGLE-2009-BLG-023/MOA-2009-BLG-028. Even though the event was 
observed solely by survey experiments, we could uniquely determine the mass of the lens and 
distance to it by simultaneously measuring the Einstein radius and lens parallax. From this, 
we find that the lens system is composed of M-type dwarfs with masses (0.50±0.07) Mq and 
(0.15±0.02) Mq located in the Galactic disk with a distance of ^ 1.8 kpc toward the Galactic 
bulge direction. The event demonstrates that physical lens parameters of binary-lens events 
can be routinely determined from future high-cadence lensing surveys and thus microlensing 
can provide a new way to study Galactic binaries. 



Key words: gravitational lensing: micro - binaries: general 



1 INTRODUCTION 

Microlensing occurs when a foreground star (lens) is closely 
aligned to a background star (source) and the light from the source 
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is refracted by the gravity of the lens. The phenomenon causes 
splits and distortions of the background stellar image. For source 
stars in the Galaxy, the separation between the split images is 
of an order of milli-arcsec and thus cannot be directly observed. 
However, the phenomenon can be observed through the bright- 
ness change of the source star caused by the change of the rela- 



tive lens-source separation ([ Paczvriski 198 ^. Currentl y, two groups 



dOGLE: Udalski et alJl2003l : lMOA: Bond et alj|2002[) are conduct- 
ing survey observations to detect microlensing events by observing 
stars located toward the Galactic bulge direction. From these sur- 
veys, more than 500 events are being detected every year. 

For most cases of Galactic microlensing events where source 
stars are lensed by a single foreground star, the lensing light curve 
is represented by 
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Here u represents the lens-source separation normalized by the Ein- 
stein radius 6e, is the time required for the source to transit 
the Einstein radius (Einstein time scale), to is the time of the clos- 
est lens-source approach, and uq is the lens-source separation at 
that moment. Among the lensing parameters of tE, uo and to char- 
acterizing lensing light curves, only the Einstein time scale pro- 
vides information about the lens because it is related to the phys- 
ical parameters of the lens mass M, relative lens-source parallax 
TTrci = AU(Dl ^ ^ ^s^)' proper motion /i by 



tE = — ; 6(e = {kMtv.^iY^^, (2) 

where k — 4G/(c^AU) and Dl and Ds represent the distances to 
the lens and source star, respectively. However, the time scale re- 
sults from the combination of the underlying physical lens param- 
eters. As a result, it is difficult to uniquely determine the individual 
lens parameters from the time scale alone. 

The degeneracy of the lens parameters can be partly lifted by 
measuring either a proper motion or a lens parallax, and can be 
completely broken by measuring both. The proper motion is re- 
lated to the Einstein radius and the time scale by /i — O^/tE 
and thus measuring /i is equivalent to measuring 6e- Einstein 
radii are generally measured from the deviation of the light curve 
from that of a point-source event caused by the finite-source ef- 



fect I Nemiroff & Wickramasinghelll994 Iwitt & Maolll994 iGouldl 



1 1994 ). The microlens parallax is defined by the ratio of the Earth's 
orbit to the Einstein radius projected on the observer plane, te, i.e., 



AU 



(3) 



Lens parallaxes are generally measured by analyzing the deviation 
of the light curve caused by the change of the observer's position 
over the course of the event due to the orbi t al motion of the E arth 
around the Sun ( lRefsdallll966l : lGouldlT99l : ISmith et al.ll2003l) . As 
a result, parallaxes are usually measured for events with long time 
scales that are comparable to a significant portion of the Earth's 
orbital period, i.e. 1 yr. If both Einstein radius and lens parallax 
are measured, the mass and the distance to the lens are uniquely 
determined by 



M = 



Be 



AU 



(4) 



KTTb TTe^E + TTS 

respectively. Here tts — AU/Dg represents the parallax of the 
source star. 

A fraction of lensing events are produced by lenses com- 



posed of binary masses jMao & Paczviisldlll99lh . These binary- 
lens events provide a good chance to determine the physical pa- 
rameters of lenses. The first reason for this is that most binary-lens 
events are detected through the channel of caustic crossings. The 
magnification gradient in the region around the caustic is so steep 
that finite-source effect is always manifested in the light curves 
of caustic-crossing events. If the caustic-crossing part of the light 
curve is resolved, then, it is possible to measure the Einstein ra- 
dius. Second, the average mass of binary lenses is heavier than that 
of single lenses and thus the time scales of binary events tend to be 
longer than those of single-lens events. This implies that the chance 
to measure lens parallaxes is higher. 

Despite these advantages, mass measurements of binary 
lenses by measuring both 8e and tte were possible only for 
a handf ul number of events . These events i nclude EROS-2000 - 
BLG-5 jAlcock et alj 1200 ll: lAn et al\ 1200 ll : iGould et alj l2004h . 

" 20051) , OGLE-2005-BLG- 
2005"; 'Dong etalj|200 9). OGLE-2005-BLG^ 
2010b), OGLE-2006-BLG-109 iGaudietal 



OGL E-2002-BLG-069 |Kubas et alj 



071 lUdalski et al 
153 (Hwang et al. 



2008 : iBennett et al.l I2010I) . MOA-2007-BLG-192 (Bennettetal 
2008I) . MOA-2009-BLG-016 



OGLE- 



Hwang et alj 

2009-BLG-092/MOA-2009-BLG-137 tevu et alj [20Tcl) . Among 
these events, three (OGLE-2005-BLG-071, OGLE-2006-BLG- 
109, MOA-2007-BLG-192) are planetary events for which inten- 
sive follow-up observations were conducted during the events. Un- 
fortunately, follow-up observations are not being conducted for 
general binary events due to the limited observational resource. 

However, the situation is rapidly changing with the instrumen- 
tal upgrade of the survey experiments. The second phase of the 
MOA survey experiment started in 2006 by replacing its old 0.6 
m telescope with a new 1.8 m telescope equipped with a camera 
of 2.2 deg^ field of view. In 2010, the OGLE group started its 
fourth phase survey observations with a new camera of 1.4 deg'^ 
field of view that is 3.5 times wider than that of the camera used 
in the third phase experiment. With the improved photometric pre- 
cision combined with the increased monitoring cadence provided 
by the upgraded instrument, the survey experiments detect much 
more events. For example, the MOA group detected 563 events in 
2009 season compared to 56 events in 2005. In addition, it became 
possible to constrain lenses based on the data obtained from survey 
observations alone for an increasing number of events. 

In this paper, we present the result of the analysis of 
a caustic-crossing binary-lens microlensing event OGLE-2009- 
BLG-023/MOA-2009-BLG-028. This event demonstrates that 
physical parameters of binary lenses can be constrained based on 
the data obtained from survey observations. We provide the physi- 
cal parameters of the lens system. We also discuss the use of future 
microlensing surveys in the studies of Galactic binaries. 



2 OBSERVATION 

Figure [T] shows the light curve of the event OGLE-2009-BLG- 
023/MOA-2009-BLG-028. As evidenced by the two strong per- 
turbations at HJD ~ 2454905 and 2454920 and the characteris- 
tic "U"-shape trough region between them, the event is a typical 
caustic-crossing binary-lens event. The brightening of the source 
star was noticed in the early 2009 bulge season by both OGLE and 
MOA survey experiments using the 1.3 m Warsaw telescope of Las 
Campanas Observatory in Chile and the 1.8 m of Mt. John Obser- 
vatory in New Zealand, respectively. The perturbation produced by 
the first caustic crossing was detected on March 12 by both survey 
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Figure 1. Light curve of the microlensing event OGLE-20()9-BLG-023/MOA-2009-BLG-028. Also presented are the model curves for the best-fit solutions 
with and without the parallax effect. We note that the MOA data except the caustic-induced perturbation part (2454900 < HJD < 2454920) are binned for 
a clear view. The zoomed view of the caustic-induced part of the light curve is presented in Fig. 4. 



experimeiits. The time gap during 2454780 <; HJD <, 2454860 
corresponds to the period during which the bulge could not be seen. 
We note that OGLE observation was stopped on HJD ~ 2454960 
for the upgrade of the camera. In addition, no follow-up observa- 
tions of the event was conducted. Nevertheless, the light curve of 
the event was well covered during the whole 2009 bulge season. 

For the analysis of the light curve, we use 1058 /-band and 15 
l^-band OGLE images taken during 2453417 < HJD < 2454955 
and 2453470 < HJD < 2454951, respectively. The^MOA data 
is composed of 1298 i?-band images taken during 2454495 <: 
HJD < 2455134. The photometry was processed by the individual 
groups using their own software. 



3 MODELING 

Light curves of binary-lensing events result from a com- 
plex phenomenology and thus exhibit an astonishing diversity 
(ISchneider & WeissHl986h . As a result, modeling light curves is a 
difficult task. One important difficulty arises due to the large num- 
ber of parameters to be included in modeling. These parameters are 
needed to describe various features of the light curve. To describe 
light curves of standard single-lens events, a set of three parameters 
of tE, tit, and no are needed. To describe the deviation caused by the 



lens binarity, an additional set of three binary parameters is needed. 
These binary parameters include the mass ratio between the lens 
components, q, the projected binary separation in units of the Ein- 
stein radius, s, and the angle of the source trajectory with respect to 
the binary axis, a. Since OGLE-2009-BLG-023/MOA-2009-BLG- 
028 is a caustic-crossing binary event, an additional parameter of 
the normalized source radius, p* = S*/Se, is needed to account 
for the finite-source effect. Here 0* represents the angular radius of 
the source star. In addition, the event lasted throughout the whole 
2009 bulge season and thus it is needed to check the possibility of 
deviations induced by the parallax effect. To incorporate the par- 
allax effect, it is required to include two parallax parameters vTE.iv 
and tte.e, which are the components of the microlens-parallax vec- 
tor tte projected on the sky in the north and east celestial coordi- 
nates, respectively, where the direction of the parallax vector is that 
of the lens-source relative motion in the frame of the Earth at the 
peak of the event. Due to the shear size of the parameter space, 
brute-force searches for solutions are very difficult and extremely 
time-consuming. 

The second important difficulty in modeling binary-lensing 
light curves is caused by the complexity of surface. This com- 
plexity implies that even if a solution that apparently describes the 
observed light curve is found, it is difficult to be sure that all possi- 
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ble minima have been investigated jPominiljl 1 999al lbl) . As a result, 
a simple downhill approach to search for solutions cannot be used. 

For efficient modeling but avoiding the difficulties mentioned 
above, we use a hybrid approach of parameter searches. In this ap- 
proach, grid searches are conducted over the space of a subset of 
parameters and the remaining parameters are searched for by let- 
ting them vary so that they result in minimum at each set of the 
grid parameters. We choose s, q, and a as the grid parameters be- 
cause they are related to the light curve features in a complex way 
such that a small change in the values of the parameters can lead 
to dramatic changes in the resulting light curve. On the other hand, 
the other parameters are more directly related to the light curve 
features and thus they can be searched for by using a downhill ap- 
proach. For the minimization in the downhill approach, we use 
a Markov Chain Monte Carlo method. Once the minima of the 
individual grid points are determined, the best-fit model is obtained 
by comparing their values. We investigate the degeneracy of so- 
lutions by probing local minima that appear in the space of the grid 
parameters. 

The last major difficulty in binary-lensing modeling arises due 
to the fact that the modeling requires large computations. Finding 
the best-fit parameters and evaluating their uncertainties requires 
generating a large number of trial model light curves. The prob- 
lem is that most binary-lensing events exhibit deviations induced by 
the finite-source effect and calculating finite-source magnifications 
requires intensive computations. Formally, the magnification of a 
finite source can be calculated by integrating the point-source mag- 
nifications of the elements of the the source star over its surface. 
However, this approach is impractical due to the divergent nature of 
the magnification near the caustic and the large amount of compu- 
tation time required for precision calculation. Although there exist 
semi-analytic approximations ( Schneider & Weiss 1986), it is not 
precise enough to describe the caustic-crossing feature of lensing 
light curves. Therefore, an efficient method of finite-magnification 
calculation is important for binary-lensing modeling. 

For the finite-source magnification calculation, we use a mod- 
ified version of the ray-shooting method. In the usual ray-shooting 
method, one shoots a large number of uniformly spaced rays in 
the image plane, and determine which ones land on the surface of 
the source using the lens equation. Then, the magnification cor- 
responding to a source position is calculated as the ratio between 
the number density of rays on the image plane to that of rays o n 
the source surface jKavser et al]|l986l : IWambsganss et al.lll99d) . 
Based on this basic scheme, we minimize the computation time by 
using the following methods. 

(i) Magnification map making 

Instead of calculating the magnifications corresponding to the in- 
dividual source positions, we construct ma gnification maps of the 
region encompassing perturbation regions Jpong et aklliooeh . The 
main advantage of this method is that once a map for a fixed (s,q) 
parameter set is constructed, one can reuse it for the calculations 
of many light curves resulting from different combinations of other 
parameters instead of re-shooting rays all over again. 

(ii) Minimization of ray-shooting region 

We further reduce the calculation time by minimizing the area of 
the ray-shooting region. We set the regions of ray shooting in the 
image plane only for which rays arrive in the region around caustics 
on the source plane where the finite-source is important. For exam- 
ple, if the perturbation occurs near the peak of a high-magnification 
map, the perturbation is localized in the central region around the 
binary lens components, where a small caustic is located. This re- 
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Figure 2. Geometry of the lens system for the best-fit parallax model. 
The two dots represent the locations of the binary lens components and the 
closed figure composed of concave curves represents the caustic. The co- 
ordinates (5, Tj) are centered at the primary lens and the abscissa is aligned 
with the binary axis. The line with an aiTow represents the source trajec- 
tory. Note that the source trajectory is curved due to the parallax effect. 
All lengths are normalized by the Einstein radius corresponding to the total 
mass of the binary. 

gion in the source plane corr esponds to the annulus around the Ein- 
stein ring in the image plane jOriest & Safazadehll998l:lDong eTai] 
I2OO6). Then, by shooting rays only in the localized region in the im- 
age plane, we minimize the number of rays needed for finite-source 
magnification calculations and thus reduce the calculation time. In 
Appendix A, we describe how the region of rayshooting is set. 

(iii) Semi-analytic approximation 
Finally, we further speed up computations by limiting numer- 
ical computation of finite-source magnifications only when the 
source is located very close to the caustic and using sim- 
ple semi-analytic hexadec apole approximations jGouldl Ia. 20081 : 
IPeicha & Hevrovs^l2009h in the vicinity of the caustic. For this, 
we divide the source plane into small grids and then register the 
individual grids by different levels depending on the distance from 
the caustic. Based on the levels, we apply different levels of ap- 
proximations for the magnification calculations. 



4 BEST-FIT MODEL 

From modeling, it is found that OGLE-2009-BLG-023/MOA- 
2009-BLG-028 is produced by the crossings of a Galactic bulge 
F-type main-sequence source star over the caustic produced by a 
disk binary lens composed of M-type main-sequence stars. The 
determined values of the normalized star-planet separation and 
planet/star mass ratio are 

s= 1.65 ±0.02, g = 0.30 ±0.01, (5) 

respectively. In Table [T] we present the lensing parameters deter- 
mined from modeling. The model light curve is presented in Figure 
[T] Figure |2] shows the geometry of the lens system under the best- 
fit model, where the two dots represent the locations of the lens 
components, the closed figure composed of concave curves repre- 
sents the caustic produced by the binary lens, and the line with an 
arrow represents the source trajectory. We note that the trajectory 
is curved due to the parallax effect. We find no serious local min- 
ima except the one caused by the mirror-image symmetry between 
the source trajectories with the impact parameters and orientation 
angles of the source trajectory of (uo,a) and (— ito, —a). We find 
that the solution with < is preferred over the solution with 
Mo > by Ax^ = 1.0, implying that it is difficult to distinguish 
the two solutions. 
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Table 1. Fit Parameters. HJD' = HJD - 2450000. 





no 


parallax 


parallax 




parallax 


(uo > 0) 


(uo < 0) 


X /dof 


4067.18/2365 


2371.12/2363 


2370.08/2363 


s 


1.860±0.015 


1.705±0.020 


1.645±0.022 


1 


0.260±0.004 


0.304±0.008 


0.295±0.012 


a (deg) 


106.87±0.15 


102.50±0.29 


-102.34±0.40 


to (HJD') 


4915.48±0.06 


4913.72±0.10 


4913.85±0.19 


MO 


0.094±0.003 


0.116±0.004 


-0.132±0.006 


t-E (days) 


153.556±2.622 


100.87±2.315 


92.49±2.61 


P* 


(2.9 ±0.3) X 10""' 


(8.4 ±0.6) X 10-* 


(9.5 ±0.7) X 10"* 






0.081±0.0I5 


-0.120±0.033 






0.231 ±0.008 


0.252±0.111 
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Figure 3. Upper panel: Contour of x in parameter space of 
(tte.Ei TTE.iv), which are the components of the lens-parallax vector tte 
projected in the sky in the north and east celestial coordinates, respectively. 
Lower panel: Contours of in the parameter space of the nonnalized 
source radius and the Einstein time scale t e . 



4.1 Lens Parallax 

Light curves of long time-scale events are susceptible to the 
parallax effect. The time scale of the event OGLE-2009-BLG- 
023/MOA-2009-BLG-028 is ~ 100 days, which comprises 
nearly 1/3 of the orbital period of the Earth. We, therefore, search 
for solutions considering the parallax effect. From this, we find that 
the parallax model provides a significantly better fit with Ax^ = 
1697. The determined values of the parallax parameters are 
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Figure 4. Enlargement of the caustic-crossing part of the light curve. The 
inset shows the light curve during the caustic entrance of the source star. 
For comparison, we draw light curves for different values of the normalized 
source radius p» . 



tte.jv 



-0.12 ±0.03, 



TTE.B = 0.25 ±0.01. 



(6) 



The upper panel of Figure |3]shows the scatter plot of in the pa- 
rameter space of (tte.b, vrE,jv). The contours are elongated along 
the 7rE,jv axis because the apparent motion of the Sun at to pro- 
jected onto the sky is perpendicular to the tte.jv axis. 

It is known that the orbital motion of a binary source can give 
rise to distortions of lensing light curves ('xalla rap' effect) that ar e 
similar to those induced by the parallax effect dSmith et al.ll2003h . 
We check this possibility by conducting xallarap modeling under 
the assumption that the binary source is in a circular orbit. From 
this, we find that the improvement of the fit of the xallarap model 



from that of the parallax model is merely 



2.4, which 



is much smaller than the improvement by the parallax effect of 
~ 1700. In addition, the minimum occurs at an orbital 
period of P ~ 1 yr, which corresponds to the orbital period of the 
Earth around the Sun. Both facts support the parallax inteipretation 
of the light curve deviation. 
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Table 2. Physical Parameters 





parallax (uq > 0) 


parallax {uq < 0) 


(mas) 


1.583±0.191 


1.475±0.187 


fi (mas yr~ ^ ) 


5.727±0.692 


5.822±0.738 


Dl (kpc) 


1.933±0.203 


1.845±0.214 


Ml (Mq) 


0.608±0.077 


0.500±0.071 


M2 (Mq) 


0.185±0.023 


0.148±0.021 



4.2 Einstein Radius 

The Einstein radius is measured from the normalized source radius 
combined with the information about the angular source radius 
0^, by 



(7) 



The normalized source radius is measured with a moderate uncer- 
tainty from the analysis of the light curve during the caustic cross- 
ings. In the lower panel of Figure |4] we present the enlargement 
of the caustic-crossing part of the light curve. From the figure, it is 
found that a single OGLE data point taken during the time when 
the source is on the fold caustic (at HJD ^ 4903.78) and multi- 
ple MOA data points taken at the time when the source is about 
to leave the caustic (during 4904.05 <, HJD ^ 4904.20) provide 
constraints on p*. The measured value of the normalized source 
radius is 



p* = (9.5 ±0.7) X 10" 



(8) 



Considering that the normalized source radius of a main-sequence 
source star is p* ~ (O)IO^^ for typical lensing events caused by 
low-mass lens stars located half way between the source and ob- 
server, the measured value of p* is substantially small. This means 
that the Einstein radius is big, suggesting that either the lens is 
heavy or it is located close to the observer. The lower panel of Fig- 
ure [3] shows the distribution in the parameter space of p* and 

The angular source radius is determined from the inform ation 
of the de-reddened color of the source star dYoo et al.ll2004l) . We 
determine the color by using the centroid of clump giant stars in 
the color-m agnitude diagram as a reference position dStanek et al.l 
Il994ll997l) under the assumption that the source star and clump 
giants experience the same amount of extinction^ The offset be- 
tween the source and clump centroid is measured in the instru- 
mental color-magnitude diagram that is constructed by using the 
OGLE V and / band images taken toward the bulge field where 
the source is located (Figure |5). Then, the angular source size is 
determined by first transfor ming from (V — I)o to {V — K)o using 
the color-color relation of iBessel & Bretj d 1998b and then apply- 
ing the relation between {V — K)o and the angular stellar radius 



^ For source stars in the bulge, the uncertainty caused by this assumption 
would be small. This is because the extinction toward the Galactic bulge 
field is caused by the dust distributed in the disk and thus even if there 
exist some differences between the distances to the clump centroid and the 
source star, the amount of extinction would be nearly the same. The source 
star can be located in the disk but the probability is low. The validity of the 
source location in the bulge can also be checked by its location on the color- 
magnitude diagram. From the diagram presented in Figure |5] it is found 
that the source is located in the region where bulge main-sequence stars are 
densely populated, suggesting that the source is located in the bulge. 
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Figure 5. Position of the source star in the instrumental color-magnitude 
diagram constructed based on OGLE data. 



of iKervella et al. I l l2004l) . For the best-fit model, we find that the 
de-reddened magnitude and color of the source star are 7o = 16.7 
and {V — J)o = 0.67, respectively, implying that the source is an 
F-type main-sequence star with an angular radius of 

£»* = (1.40 ± 0.15) fias. (9) 

Here we adopt an average distance to clump giants toward the field 
of 7.7 kpc estimated ba sed on the Galactic mass distribution model 
of iHan & Goui3 ilOOj ). The uncertainty of 6** is determined from 
the combination of the uncertainty of the distance to the source and 
an additional 7% intrinsic error in the conversion process from the 
measured color to the source radius. 

With the measured values of the normalized and angular 
source radii, the Einstein radius is determined as 

= (1.48 ±0.19) mas. (10) 

This corresponds to the relative lens-source proper motion of 



(5.82 ±0.74) masyr'\ 



(11) 



4.3 Pliysical Parameters 

With the measured Einstein radius and parallax, the mass and dis- 
tance to the lens are determined from the relations in equation l|4]l. 
For the best-fit model, these values are 



M = (0.65 ± 0.09) Mq 
and 

Dl = (1.8 ±0.2) kpc, 



(12) 



(13) 



respectively. With the determined mass ratio, it is found that the 
masses of the individual lens components are Mi — (0.50 ± 
0.07) Me and M2 = (0.15 ± 0.02) Mq, respectively. Therefore, 
the lens is composed of an early and a late M-type main-sequence 
stars located in the Galactic disk. The model with ito > yields 
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the distance and mass of the lens slightly bigger than those of the 
model with uq < (see Table|2ll. 



5 DISCUSSION AND CONCLUSION 

We analyzed the light curve of the long time-scale, caustic- 
crossing, binai-y-lens event OGLE-2009-BLG-023/MOA-2009- 
BLG-028. Despite that the event was observed solely by survey 
experiments, we could uniquely determine the mass of the lens and 
distance to it by simultaneously measuring the Einstein radius and 
lens parallax. It was turned out that the event was produced by the 
crossings of a Galactic bulge F-type main-sequence source star over 
the caustic produced by a disk binary lens composed of M-type 
main-sequence stars. The event demonstrates that the physical pa- 
rameters of binary lenses can be uniquely determined from data 
obtained by survey observations. 

Microlensing can potentially probe the distributions of bi- 
naries as functions of mass ratio and separation, that can pro- 
vide important observational constraints on theories of star for- 
mation ((Gould 2001). Especially, microlensing is sensitive to low- 
mass companions that are difficult to be detected by other meth- 
ods and thus it is possible to make complete distributions down 
to lower mass limit of binary companions. Due to this impor- 
tance, there have been several systematic searches fo r binary lenses 
dAlcock et al. 2000; Jaroszvhski 2002; Jaroszvnsk i. et al.l |2004| . 



20051 l2006l ; 1skowron et al.ll200 7). However, the samples of bina- 



ries acquired from previous surveys were not adequate enough to 
strongly constrain the binary distributions. One important reason 
for this is the difficulty in estimating the detection efficiency of 
binary-lens events. Most binary-lens events are detected through 
the channel of caustic-crossing events where caustic crossings were 
accidently discovered by the sudden rise of the source star flux. 
Due to this haphazard nature of binary-lens events, it is difficult to 
estimate the detection efficiency that is essential for the statistical 
studies of binaries. In addition, the physical quantities of lenses for 
most of binary-lens events could not be determined, making de- 
tailed studies of binaries difficult. 

However, the situation will be different with the advent of 
new-generation experiments. As mentioned, the recent upgrades of 
the OGLE and MOA experiments already significantly increased 
the observational cadence of the surveys. In addition to these exper- 
iments, there is a planned experiment that can increase the cadence 
even higher. Korea Microlensing Telescope Network (KMTNet) 
is an approved project that will employ three telescopes, each of 
which will have a 1 .6 m aperture and 4 deg^ field of view. They will 
be located in three different continents of South America (Chile), 
Africa (South Africa), and Australia for continuous observations of 
microlensing events. The expected cadence of the experiment is 6 
hr~^. Considering that the typical time scale of caustic crossings 
is several hours for events involved with main-sequence stars and 
extends to > 10 hrs for events associated with giant source stars, 
the cadence of the survey is high enough to resolve caustic cross- 
ings of most events, enabling measurements of Einstein radii. For 
a significant fraction of these events, it will be possible to addi- 
tionally measure lens parallaxes, enabling to completely determine 
the physical parameters of binary lenses. Furthermore, the uniform 
coverage of events will make it easier to estimate the detection effi- 
ciency of binary-lens events, enabling statistical analysis of binary 
distributions possible. Therefore, future lens surveys can provide a 
new way to study Galactic binary stars. 
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APPENDIX A: REGION OF RAYSHOOTING 

We determine the range of rayshootinh as follows. For sources 
close to a fold caustic, the point-source magnification is represented 
by 

A^^ Ao + {d/u^r^'''(d{d), 

where d is the perpendicular distance to the caustic with d > 
for sources interior to the caustic, u±_ represents the characteristic 
strength of the local caustic, 0(a;) is the Heaviside step function, 
and Ap represents the total magnification of the s lowly varying im- 
ages dSchneider & Weiss ' 198a lOaudi & F etters 2002). Then, un- 
der the approximation that the magnification variation of the slowly 
varying images is negligible, the point-source approximation can 
be used for magnifications in the region outside the caustic with 
|d| > p*. Inside the caustic, the finite-source magnification, Ai, is 
expressed as 

/P* / 2 _ 2 \ 1/2 
(^rf j * = d-''\f{d/p.). 

Then, the fractional deviation from the point-source magnification 
is 

Ai - Ar, 



implying that the fraction deviation depends only on the ratio d/p*. 
We find that the deviation are 2%, 1%, and 0.5% for the distances 
from the caustic of d = 2.3p*, 3.1p*, and , 4.3p*, respectively, 
implying that the deviation decreases rapidly with the increase of 
d/p*. To be conservative, we set the range of rayshooting as d < 
2p* and d < 5p* inside and outside the caustic, respectively. 



